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ABSTRACT 

The rapid, irreversible change of the photospheric magnetic field has been recognized as an important element 
of the solar flare process. This Letter reports such a rapid change of magnetic fields during the 201 1 February 
13 M6.6 flare in NOAA AR 11158 that we found from the vector magnetograms of the Helioseismic and 
Magnetic Imager with 12-min cadence. High-resolution magnetograms of Hinode that are available at ~-5.5, 
-1.5, 1.5, and 4 hrs relative to the flare maximum are used to reconstruct three-dimensional coronal magnetic 
field under the nonlinear force-free field (NLFFF) assumption. UV and hard X-ray images are also used to 
illuminate the magnetic field evolution and energy release. The rapid change is mainly detected by HMI in a 
compact region lying in the center of the magnetic sigmoid, where the mean horizontal field strength exhibited 
a significant increase by 28%. The region lies between the initial strong UV and hard X-ray sources in the 
chromosphere, which are cospatial with the central feet of the sigmoid according to the NLFFF model. The 
NLFFF model further shows that strong coronal currents are concentrated immediately above the region, and 
that more intriguingly, the coronal current system underwent an apparent downward collapse after the sigmoid 
eruption. These results are discussed in favor of both the tether-cutting reconnection producing the flare and 
the ensuing implosion of the coronal field resulting from the energy release. 

Subject headings: Sun: activity — Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: X-rays, gamma 
rays — Sun: magnetic topology — Sun: surface magnetism 



1. INTRODUCTION 

^ \ It has been known that the long-term evolution of photo- 
OO spheric magnetic field (PMF) driven by new flux emergence 
On and surface flows plays important roles in building up free 
energy in the corona, and that this free magnetic energy pow- 
. ers flares and coronal mass ejections (CMEs; Priest & Forbei 
(S| ■ [2002). On the other hand, a short-term variation of the PMF 
^-H ' associated with flares has not been considered because PMFs 
^-H are strongly line-tied to the dense high-/3 photosphere and 
^-H thus are thought unlikely to be altered by any flare-related 
disturbances created in the tenuous low-/3 corona. Only re- 
. ^ cently, a back reaction of the coronal magnetic field (CMF) 
^ on the PMF during the reconfiguration of the CMF has been 
^ seriously considered from the theoretical point of view. The 
pi idea is that the CMF should contract inward, as the mag- 
netic energy of the CMF decreases after flares and/or CMEs 
(called "implosion"; H udsonll200 0). This may make the PMF 
be oriented more horizontally, resulting i n a Lorentz force 
acting dow nwardly on the solar surface (iFisher et al.l 120101: 
[Hudson et al.> 2008). In fact, a few theoretical flare mod- 
els predict the appearance of more horizontal PMFs after 
flares/CMEs in the form of ne wly fornied, lo ^ y-lying field 
lines close to the surface (e.g., lMelrosd[l997l: iMoore et aH 

Hooil). 

The phenomenon of coronal implosion has been re- 
ported in several recent o bservational studies (iLiu et al .1120091: 
iLiu & Wangl l2009L 12010 ). The related phenomenon of the 
rapid change of the PMF to incline more horizontally after 
flares had also been repo rted almost two decades ago ( WangI 
1992: Wang et al.l 11994 ). Namely, the transverse field near 
the magnetic polarity inversion line (PIL) at the core flar- 



ing region is enhanced rapidly and irreversibly, which is of- 
ten accompanied by an increase of magnetic shear. Such 
a rapid change of vector PMF change closely associated 
with flares/CMEs has been consistently found later on us- 
ing ground-b ased jS chmieder et al."1994'; "Wang etal." "20021 
20p4b, 2005] ILiu et al. 2005; Wang et al. 2007; Wang & Liu 
20101) and Hinode observations ( Jing et al. 2008; Li et al.| 
20091) . Indirect evidence includes the unbalanced fl ux evo - 



lution of the line-of- sight magnetic field (I Wang & Lir 
the variation of the sunspo t white-light struc t ure in flaring re- 
gions iWans et al. 2004a; iDeng et al.l 120051: iLiu et al. 1 120051: 
IChen et al.» >2007). and the change in the pattern of the penum- 
bral Ever shed flow (Deng et al. 2011). Additional support 
also comes from the three-dimension al (3D) M HD simula- 
tions of an emerging twisted flux tube CFanl2010|) . the eruption 
of whi ch induces sim ilar impact on the low-altitude magnetic 
fields ( Li et al.ll20rih . 

While the evolution of vector PMFs closely associated with 
flares were mostly studied using ground-based observations 
subjec t to seeing variations, Wang et al. (2011) and Sun et alj 
(1201 1|) use d data from the H elioseismic and Magnetic Im- 
ager fHMI: ISchou et al.|[201 1|) on board the recently launched 
Solar Dynamics Observatory (SDO) to find a rapid enhance- 
ment of transverse field associated with the 201 1 February 15 
X2.2 flare. Since solar activity is now increasing, this is an op- 
portune moment to advance the study of the flare-related evo- 
lution of the PMF using the state-of-the-art, seeing-free pho- 
tospheric vector magnetograms acquired by space missions. 

In this Letter, we investigate the PMF evolution associated 
with an M6.6 flare on 201 1 February 13 that was well covered 
by both HMI and Hinode (see Section O, in which the trans- 
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Figure 1. By maps are superimposed with arrows (for clarity, those in negative and positive fields are coded in different colors) representing horizontal plasma 
flows (a) and magnetic field vectors (b-d). The NLFFF lines in (e) stem from the four flare kernels FP1-FP4 in 1700 A (f) at the flare onset. The white contours 
in (e) outline the double J-shaped flare ribbons in Ca ii H. Contours (with levels of 38%, 48%, 58%, 68%, 78%, 88% of the maximum flux) in (f) denote RHESSI 
CLEAN images reconstructed with detectors 2-8 (except that 17:32:08 UT image is made without the detector 2 due to the relatively low counts at this time). 
The white box overplotted on all panels denotes the region R under study with enhanced Bh after the flare. 
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Figure 2. Temporal evolution of (Bh) of the region R in Figure[T] in comparison with the flare light curves of soft and hard X-rays. The mean error of of the 
region R at each time instance measured by HMI is plotted as error bars. The error estimates for Hinode data are currently not usable^. 



verse field enhancement at the flaring PIL is unequivocaUy 
detected. We wiU discuss the rapid change of the 3D CMF as- 
sociated with that of the PMF using nonhnear force-free field 
(NLFFF) extrapolations from the active region. 

2. OBSERVATIONS AND DATA REDUCTION 
The HMI instrument obtains filtergrams in six polarization 
states at six wavelengths along the Fe I 6173 A spectral line 
to compute Stokes parameters IQUV. The inverted and disam- 
biguated vector magnetic field data for the NOAA AR 11158 
with a 12 minute cadence and 0.5'^ pixel size were recently 
released by the HMI tearrQ ( Hoeksema et al. 201 1). We used 
the vector magnetograms remapped to heliographic coordi- 
nates to (1) investigate the rapid PMF evolution associated 
with the M6.6 flare on 2011 February 13, which started at 
17:28 UT, peaked at 17:38 UT, and ended at 17:47 UT in 
terms of GOES 1-8 A flux, and (2) trace the surface flows 
with the differential affine velocity estimator for vector mag- 
netograms (DAVE4VM; Schuck 2008). 

The Spect ro-Polarimeter (SP) of the Solar Optical Tele- 
scope (SOT; iTsuneta et al.l ^2008) on board Hinode obtains 
polarization profiles of two magnetically sensitive Fe lines 
at 630.15 and 630.25 nm to generate Stokes images, which 
are then reduced with the Hinode/SF calibration/inversion 
pipeline to retrieve the vector PMF using the Milne-Eddington 
gRid Linear Inversion Network (MERLINfl Hinode/SF 
scanned the nearly entire isolated, near disk-center (about 
S20°, E04°) AR 11158 at four time bins (12:00-12:32, 
16:00-16:32, 19:00-19:32, and 21:30-22:02 UT) close to 
the M6.6 flare. We further processed the resulted vector 
magnetograms with a resolution of 0'^32 pixel"^ to (1) re- 
solve the 180° azimuthal ambiguity using the AZAM code 
(Leka et al. 2009) based on the "minimum energy" approach 
( Metcalf et al. 2006), and (2) remove the projection effects by 
transforming the observ ed fields to heliographic coordinates 
(iGary & Hagyardll990|) . Taking advantage of the high spatial 

^|http : // jsoc . stanford.edu/jsocwiki/VectorPaper 

^ http : / /sot . 1 ms al . co m/data/sot /level2dd/ 



sot sp_level2_de sc ript io n .html 



resolution and high polarization accuracy of Hinode, we ex- 
amine the CMF evolution by constructing NLFFF models us- 
ing the "weighted optimization" method ( Wiegelmann 2(j04|) 
after preprocessing the photospheric boundary t o best suit 
the force-free condition (IWiegelmann et al.l 12006 ^. The cal- 
culation was performed using 2x2 rebinned magnetograms 
within a box of 236 x 256 x 256 uniform grid points, which 
corresponds to ~110 x 120 x 120 Mm^ of balanced mag- 
netic fluxes (with a ratio of opposite surface fluxes of ~99%- 
102%). 

The temporal and spatial relationship between the PMF 
change and flare energy release can provide important clues 
concerning the eruption mechanism. The evolution of the 
flare hard X-ray (HXR) emission was entirely registered by 
the Reuve n Ramaty Hi^h E nergy Solar Spectro scopic Imager 
{RHESSI\^m^m^^- CLEAN images (iHurford et al.l 
2002) in the nonthermal energy range (25-100 keV) show- 
ing the flare footpoints were reconstructed using the front seg- 
ments of detectors 2-8 with '^48 s integration time throughout 
the event. To provide the observational context from the low 
chromosphere to the coronal, we also used 1700 A continuum 
(5000 K) and 94 A (Fe xviii; 6 MK) images taken by the 
Atmospheric Imaging Assembly (AIA: Lemen et al.ll201 Ih on 
board SDO, and Ca II H images taken by Hinode/SOT. 

3. FLARE-RELATED PMF CHANGES 

The M6.6 flare is initiated at the center of the active region, 
where opposite magnetic flux concentrations undergo an over- 
all counterclockwise, rotation-like motion as clearly tracked 
by DAVE4VM (FigurellJa)), resulting in highly sheared fields 
along the PIL (Figures \V[h) and (c)). A compact region R 
at the central PIL (the boxed region in Figure [B is readily 
identified in the HMI vector magnetograms to have a pro- 
nounced, ~28% enhancement of the mean horizontal field 
strength from (5/,) c^l 160 G at 17:24 UT immediately before 
the flare to {Bh) c^USO G at 18:00 UT right after the flare 
(see Figure [21) in ~30 minutes, with a rate of ~0.15 G s~^. 
The change-over time is cotemporal with the rapid rising of 
soft X-ray flux and peaking of HXR emissions. Standing out 
over the long-term evolution of (Bh), which increases much 
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(a) Magnetic shear 
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Figure 3. (S), S, (cp), and (Jh) of the region R in Figure [T] as a function of 
altitude with a step size of ~0.46 Mm, at four times of Hinode scans. 

more slowly with time (at ~8 x 10~^ G s~^), this change oc- 
curs in a stepwise fashion, coincident with the occurrence of 
the M6.6 flare, and the enhanced {Bh) survives the flare and 
persists pernianentl y. Consistent with our previous results 
(IWang & Liull2QlQl and references therein), this rapid PMF 
change also leads to a more horizontal state of the PMF in the 
region R after the flare. 

Similar to the HMI result, the {Bh) within the region R is 
observed by Hinode to increase from c^l090 G (12 UT) and 
1180 G (16 UT) before the flare to ::^1600 G (19 UT) and 
1660 G (21:30 UT) after the flare (cf. Figures [Hb) and (d), 
and (g) and (h); also see Figure [2]). It is particularly inter- 
esting that there are four compact flare kernels FP1-FP4 in 
1700 A at the flare onset, and that the region R is sandwiched 
betwee n the much stron ger FP2 and FP3 (Figure [TJf))- Fol- 
lowing iLmetal] (l2010a l). we use the locations of the flare ker- 
nels together with a preflare NLFFF model at 16:00 UT based 
on the Hinode data to shed light on the configuration of the 



flaring magnetic fields. The extrapolated field lines connect- 
ing among FP1-FP4 as plotted in Figure [TJe) apparently form 
a S-shaped sigmoidal structure with right-bearing loops con- 
necting FP1-FP2 and FP3-FP4, which is similar to the over- 
all shape of coronal loops in EUVs. This is consistent with 
the calculation that the magnetic he licity in this active region 
has a positive sign (iLiu et al.ll2011h . It is furthermore notice- 
able that (1) the early HXR emission before -17:32:30 UT 
is imaged as footpoint-like sources cospatial with FPl and 
FP2 (Figure [nO), and (2) ribbon-like HXR emissions mim- 
icking flaring loops that connect among FP1-FP4 dominate 
later on, when flare kernels evolve to double J-shaped ribbons 
(the white contours in Figure [TJe)). These lead us to a picture 
in which the flare could be triggered by the te ther-cutting re- 
connection (iMoore et al.ll 19951 1200 11: IXu et al.ll2010l) between 
the two sets of sigmoidal loops FP1-FP2 and FP3-FP4. The 
reconnected large-scale fields FP1-FP4 are seen to erupt out- 
ward to become the CME, and the newly formed smaller loops 
FP2-FP3 lying close to the surface could then account for 
the detected enhancement of Bh at the region R. We also note 
that (1) the aforementioned tether-cutting reconnection could 
be driven by the converging flows from FP3 toward FP2 (see 
FigurelUa); Moore et al. 2001; Liu et al. 2010b), and (2) the 
region R at the PIL lies between the flare kernels/ribbons (see 
Figures [He) and (f)), hence the observed Bh enhancement is 
unlikely to be affected b y flare emissions (iPatterson & ZirinI 
[19811: [Qnr&G^l2003h . 

4. CHANGE OF CMF TOPOLOGY INFERRED FROM NLFFF MODEL 

To quantify the flare-induced change of the CMF in 3D, we 
first define several characteristic quantities of magnetic field. 
The nonpotentiality can be evaluated using magnetic shear S, 
which is com puted here as th e product of field strength and 
shear angle (iWang et al.l[T99l [20061: IJing et al.ll2008D : 



S = BO 



(1) 



where B = \B\, = cos~^(B - Bp)/(BBp), and the subscript p 
represents th e potential field deri ved using the Green's func- 
tion method (iMetcalf et al .112008]) . The weighted shear angle 
over a region of interest with n pixels is then 



S = 



where / = 1,2, . 



(2) 



In this study, the magnetic inclination angle (f is defined as 
the azimuthal difference relative to the horizontal plane: 



if = tan 



(52+«2)l/2 



(3) 



Following "HudsonetaD (200S) and lFisher et'IH (l2010h . the 
change of Lorentz force acting on the surface in the vertical 
direction given a change of PMF SB can be quantified as 



downward "~ 



47r 



dA(-B^SB^ - BySBy + B^^) . (4) 



We calculate (S), S, and {(p) of the region R as a function 
of altitude, and make the comparison among the four NLFFF 
extrapolations using Hinode scans in Figures [3a)-(c). It is 
clearly shown that (1) after the flare, (S) becomes larger from 
the surface up to ~ 1 Mm but smaller above that height; mean- 
while, S decreases at nearly all heights. These suggest that 
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Figure 4. (a) A preflare Hinode By map superimposed with the NLFFF lines colored according to |/v| at footpoints. (b) A prefiare /y map superimposed with 
isosurfaces of = 0.05 A m~^ (gray), (c-f) Pre- and postflare distribution of in a vertical cross section with its bottom side denoted in (b) as the dotted line. 
The contour levels are 0.01, 0.02, 0.03, 0.05, and 0.09 A m-^. 



while the CMF is relaxed in general following the flare en- 
ergy release, the near-surface fields become more stressed due 
to a concentration of currents immediately above the region R 
(see below) with the new connectivity between FP2 and FP3 
as a result of the tether-cutting reconnection. Similarly, both 
(S) and S enhance near the surface in the event of Jin g et alJ 
(|2008), which was, however, attributed to the newly emerged 
sheared fields after the flare. (2) the whole {(p) profile shows 
increase at nearly all heights in the preflare and postflare states 
presumably due to the long-term evolution driven by new flux 
emergence. Noticeably, below ~6.5 Mm {(p) decreases across 
the flare time with larger drop nearer the surface. This clearly 
indicates that the magnetic field of and at low altitudes above 
the region R turns more inclined toward the surface, which 
agrees with the implication from the (Bh) enhancement at 



the surface and may evidence the coronal implosion scenario 
(lHudsonll20QQl) in which the CMF collapses downward after 
releasing energy. The implosion may result in a downward 
impulse due to the change of Lorentz force SF^ acting on the 
region R, which is estimated to be -2 x 10^^ dynes. 

In order to visualize the collapse of the CMF using the 
NLFFF model, we further derive the total electric current den- 
sity and its horizontal component \Jh \ = (J^-\-Jy)^^^. We note 
the followings based on the results presented in Figure IH 

First, associated with the sigmoidal fields at the flare core 
region (highly sheared field shown in the red-yellow color 
with overarching higher fields), strong currents flow along the 
sheared flaring PIL. The channel-like portion of the highest 
density //^ > 0.05 A m~^ is located directly above the region 
R (with height <1.5 Mm; Figures |4l;b)-(d)). This smaU and 
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Figure 5. Isosurfaces of Jh = 0.03 A viewed from 30° relative to the 
horizontal plane and 45° clockwise about the vertical direction. The slice is 
the same as used in Figured 



low-lying volume with the concentration of current is presum- 
ably the site of the tether-cutting reconnection ( Moore et al] 
I2OOI). It is worthwhile to point out that the most intense, ini- 
tial flare kernels FP2 and FP3 as well as the early HXR emis- 
sion centered at FP2 are not cospatial with the maximum ver- 
tical electric current density Jy but instead positioned on the 
two sides next to the above intense horizontal current channel 
(Figures [Jf) andH^b)). This suggests that it is the instability 
of the horizontal current that is responsible for electron acce l- 
eration during the early phase of the flare ( Lek a et aPI 19931) . 

Second, we cut a slice through the 3D CMF (with its sur- 
face intersection denoted in Figure Ub) as the dotted line) 
and plot the distribution of Jh in this cross section for the 
pre- (Figures HJc) and (d)) and postflare (Figures HJe) and (f)) 
states. The downward collapse of the CMF after the flare is 
clearly demonstrated by the evolution of Jh, which is most 
discernible as the contour level of 0.03 A m~^ (white bor- 
dered, dotted lines) decreases from ~3.6 to 2.3 Mm above 
the region R. Comparison between the pre- and postflare 3D 
isosurfaces of Jh (Figure [5]) also evidently indicates that the 
collapse predominantly occurs above the region R (pointed to 
by the red arrow) where the energy is expected to release as 
above mentioned. We further plot (Jh) in the region R as a 
function of altitude in Figure [3d), in which it can be seen 
that (1) (Jh) becomes stronger at nearly all heights from 12 
to 16 UT, which could be due to the long-term evolution of 
the active region (e.g.. Figures [2]), (2) from 16 to 19 UT, (Jh) 
further enhances below ~1.5 Mm with larger increase nearer 
the surface while it tends to diminish above that height, and 
(3) (Jh) does not evolve significantly from 19:00 to 21:30 UT. 
Importantly, the enhancement of {Jh) across the flare dura- 



tion is consistent with tha t of (S) at low a ltitudes, and is also 
in line with the model of iMelrosd (Il997h . in which smaller 
current-carrying loop is produced by the tether-cutting-like re- 
connection between the initial larger loops so that the current 
path moves downward closer to the surface. We tentatively 
suggest that the tether-cutting reconnection induces the CMF 
reconnection, which is accompanied by the ensuing collapse 
of the CMF due to energy release. 

5. SUMMARY AND DISCUSSION 

We have used the vector magnetograms from HMI and Hin- 
ode to analyze the change of the PMF associated with the 
2011 February 13 M6.6 flare, and interpreted the results with 
the aid of NLFFF modeling based on the Hinode data under 
the context of 3D magnetic reconnection. Major results are 
summarized as follows. 

1 . A compact region R at the central flaring PIL shows a 
rapid and permanent enhancement of {Bh) by ~320 G 
(~28% of the preflare magnitude as observed by HMI) 
closely associated with the flare. A similar trend of {Bh) 
evolution is also detected by Hinode. Overall, the PMF 
at the region R becomes more inclined to the surface 
after the flare. 

2. The magnetic fields lying immediately above R, which 
correspond to the central part of the sigmoidal region, 
are highly sheared and have strong horizontal elec- 
tric currents. The most intense flare kernels are ini- 
tially located at the two ends of this channel-like cur- 
rent system, implying that the start of electron ac- 
celeration could be closely associate d with the strong 
horizontal currents (iLeka et all 1 19931) . The configura- 
tion of the NLFFF lines is favorable for a sub sequent 
tether-cutting reconnection (iMoore et al.ll2001h . when 
disturbed by converging surface flows. This ends up 
with the short and low-lying loops, which explains the 
enhanced Bh and S at the region R. 

3. Across the flare energy-release time, the nonpotential- 
ity as represented by S increases in the height range 
from the surface up to ~1 Mm and then decreases at 
higher altitudes, which is closely correlated with the 
evolution of Jh. This is in accordance with t he recon- 
nection of current-carrying loops modeled by IMelrosd 
{1997). Remarkably, the vertical profiles of Jh as well 
as (f before and after the flare clearly demonstrate a 
downward collapse of the CMF by ~1.3 Mm toward 
the su rface , which is probably related to coronal implo- 
sion (lHudsonll2000l) . 

In summary, the change of the PMF and the corresponding 
change of the CMF observed in this study are not caused by 
the long-term development of the active region but occurred 
as rapid as the flare. The magnetic field change may involve 
two closely related physical processes: the tether-cut recon- 
nection producing the flare and the ensuing collapse of the 
CMF resulting from the energy release. 
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